Continuum percolation of overlapping discs with a distribution of radii having a 

power-law tail 
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We study continuum percolation problem of overlapping discs with a distribution of radii having a 
power-law tail; the probability that a given disc has a radius between R and R + dR is proportional 
to R~^°''^^\ where a > 2. We show that in the low-density non-percolating phase, the two-point 
function shows a power law decay with distance, even at arbitrarily low densities of the discs, uidike 
the exponential decay in the usual percolation problem. As in the problem of fluids with long-range 
interaction, we argue that in our problem, the critical exponents take their short range values for 
a > 3 — r]sr whereas they depend on a for a < 3 — r]ar where rjsr is the anomalous dimension for the 
usual percolation problem. The mean-field regime obtained in the fluid problem corresponds to the 
fully covered regime, a < 2, in the percolation problem. We propose an approximate renormalization 
scheme to determine the correlation length exponent v and the percolation threshold. We carry out 
Monte-Carlo simulations and determine the exponent v aa a. function of a. The determined values 
of V show that it is independent of the parameter a for a > 3 — rjsr and is equal to that for the 
lattice percolation problem, whereas u varies with a for 2 < a < 3 — rjsr- We also determine the 
percolation threshold of the system as a function of the parameter a. 

PACS numbers: 64.60.ah, 64.60.De, 02.50.Ey, 05.10.Ln, 05.70.Fh 



I. INTRODUCTION 

In problems like effective modeling of random media, 
the continuum models of percolation are more realistic 
than their lattice counterparts. So, much effort has been 
put into the study of such systems in the recent past. In 
two dimensions, the model systems studied involve discs, 
squares etc. of the same size or of different sizes [ll-@| 
and in three dimensions spheres, cubes etc., distributed 
randomly in space 04111 ■ An interesting sub-class of 
problems is where the basic percolating units have an 
unbounded size distribution. These are comparatively 
less studied, though a few formal results are available 
(l^ . The problem of disc percolation where discs have 
bounded sizes has been studied a lot, mainly by simu- 
lation ^ il3i, il4j] . For the single sized disc percolation, 
threshold is known to a very high degree of accuracy [11] . 
Also simulation studies have shown that the disc perco- 
lation in 2D with discs of bounded size falls in the same 
universality class as that of lattice percolation in 2D [l^ . 
For a review of continuum percolation see p^ . 

In this paper we consider continuum percolation model 
of overlapping discs in 2D where distribution of the radii 
of the discs has a power-law tail. We address questions 
like whether the power- law tail in the distribution of radii 
changes the critical behavior of the system, and how does 
the percolation threshold depend on the power of the 
power-law tail. The power-law distribution of the radii 
makes this system similar to the Ising or fluid system 
with long-range interactions. For the latter case, it is 
known that the long-range nature of the interaction does 
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affect the critical behavior of the system for slow enough 
decay of the interaction [17]. For similar results in the 
context of long-range epidemic processes, see [l^. 

The plan of this paper is as follows: In section [III we 
define the model of disc percolation precisely. In sec- 
tion mil using a rigorous lower bound on the two-point 
correlation function, we show it decays only as a power- 
law with distance for arbitrarily low coverage densities. 
We discuss the two-point function and critical exponents. 
In section IIVI we propose an approximate renormaliza- 
tion scheme to calculate the correlation length exponent 
u and the percolation threshold in such models. In sec- 
tion El we discuss results from simulation and section IVll 
contains some concluding remarks. 



II. DEFINITION OF THE MODEL 

We consider a continuum percolation model of over- 
lapping discs in two dimensions. The number density of 
discs is n, and the probability that any small area element 
dA has the center of a disc in it is ndA, independent of 
all other area elements. For each disc, we assign a radius, 
independently of other discs, from a probability distribu- 
tion Proh{R). We consider the case when Proh{R) has 
a power-law tail; the probability of radius being greater 
than R varies as i?^" for large R. For simplicity, we con- 
sider the case when radii take only discrete values RqK^ 
where j = 0,1,2,..., with probabilities (1 — p)p-' where 
p — A^°. Here i?o is the size of smallest disc, and A is a 
constant > 1. We call the disc of size RqK^ as the disc 
of type j. 

The fraction of the entire plane which is covered by at 
least one disc, called the covered area fraction /covered. 
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is given by 

/covered = 1 - CXp (1) 

where A is the areal density - mean area of the discs per 
unit area of the plane - of the discs, which is finite only for 
a > 2. For a < 2, in the thermodynamic limit all points 
of the plane are eventually covered, and f covered = 1. If 
a > 2, we have areal density, 

A = 71^i?g(l-p)/(l_pA2) (2) 

We define the percolation probability Poo as the prob- 
ability that a randomly chosen disc belongs to an infinite 
cluster of overlapping discs. One expects that there is a 
critical number density n* such that for n < n*, Poo is 
exactly zero, but Poo > 0, for n > n* . We shall call the 
phase n < n* the non-percolating phase, and the phase 
n > n* as the percolating phase. 

It is easy to show that n* < oo. We note that for 
percolation of discs where all discs have the same size 
-Ro: there is a finite critical number density n^, such that 
for n > nl, Pqo > 0. Then, for the polydisperse case, 
where all discs have radii Rq or larger, the percolation 
probability can only increase, and hence n* < n*. Also 
it has been proved that when ever we have a bounded 
distribution of radii of the discs, the critical areal density 
is g reater than that for a system with single sized discs 
. Our simulation results show that this remains valid 
for unbounded distribution of radii of the discs. 



III. THE NON-PERCOLATING PHASE 

We define two point function Prob(l 2) as the 
probability that points Pi and P2 in the plane are con- 
nected by overlapping discs. Then, by rotational invari- 
ance of the problem, Prob(l ^ 2) is only a function of 
the euclidean distance ri2 between the two points. Let 
Prob'^''(l 2) denote the probability that there is at 
least one disc that covers both Pi and P2 . Then, clearly, 

Prob(l 2) > Prob(^)(l 2). (3) 

It is straightforward to estimate Prob^^''(l 2) for our 
model. Let j be the minimum number such that radius of 
disc of type j is greater than or equal to ri2, i.e. RqA^ > 
ri2- Let S be the region of plane such that the distance 
of any point in S from Pi or P2 is less than or equal to 
i?o A-' . This region S is greater than or equal to the region 
where each point is within a distance ri2 from both Pi 
and i-2- Using elementary geometry, the area of region 
S is greater than or equal to (27r/3 — V3/4)rj2 (See Fig. 
[T]). The number density of discs with radius greater than 
or equal to RqK^ is nA~°^. Therefore, the probability 
that there is at least one such disc in the region S is 
1 — exp (— |S'|A~°-'), where \S\ is the area of region S. 
Thus we get, 

Prob(^)(l 2) > 1 - exp [~nKK-''^rl^] (4) 




Figure 1: Points 1 and 2 in the plane at a distance r from 
each other will be covered by a single disc of radius R, if the 
center of such a disc falls in the area of intersection of two 
circles with radius R and centers at 1 and 2. 



where if = 27r/3 - ^/3/4. 

Now, clearly, i?oA-' < ri2A. Hence we have A""-' > 
/ Rq" . Putting this in Eq. H we get 

Prob(^'(l w 2) > 1 - exp [-nKK-''r^^+'^] (5) 

where some constant factors have been absorbed into K. 
For large ri2, it is easy to see that this varies as r^^"' ■ 
Hence the two-point correlation function is bounded from 
below by a power-law. 

We can extend this calculation, and write the two- 
point correlation function as an expansion 

00 

Prob(l 2) = ^ Prob^") (1 ^ 2) (6) 

n=l 

where Prob'"-'(l 2) is the probability that the path of 
overlapping discs connecting points Pi and P2 requires 
n discs. The term n — 2 corresponds to a more compli- 
cated integral over two overlapping discs. But it easy to 
see that for large ri2, this also decays as r^2'^'^ ■ Assum- 
ing that similar behavior holds for higher order terms as 
well, we expect that for all non-zero densities n, the two- 
point correlation function decays as a power law even for 
arbitrarily low densities of discs. 

We note that this is consistent with the result that 
for continuum percolation in d dimensions, the diame- 
ter of the connected component containing the origin say 
{D) is divergent even for arbitrarily small number den- 
sities when {R'''^^) is divergent [11] • Here R denote the 

radii variable. In our case (D) = f ri2 ——dri2 ^ 

dri2 

J r^2'^dfi2 (where Pi is the origin) is divergent when 
a < 3, consistent with the above. 

The power-law decay of the two-point function is the 
result of the fact that for any distance r, we have discs 
of radii of the order of r. However for large values of r, 
we can imagine that there would also be a contribution 
from a large number of overlapping discs of radii much 
smaller than r connecting the two points separated by 
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the distance r, which as in the usual percolation problem 
decays exponentially with distance. Therefore it is rea- 
sonable to write the two point function in our problem 
as a sum of two parts; the first part say Gsr(?') due to 
the 'short range' connections which has an exponential 
decay with distance for large r and the second one say 
Gir{r) due to the 'long range' connections which has a 
power law decay with distance. Therefore, 



Gir) = Gsr{r) + Gir{r) 



(7) 



where 



Gir{r) ^ D{A)/r°- + higher order terms (8) 

where D{A) is assumed to go to a non-zero constant as 
A ^ A* and its dependence on A is a slowly varying one. 

Similar scenario exists whenever we have long range 
interactions in a system, such as in Ising model with 
lon g ra nge potentials or fluid with long range interactions 
[20ll2l|. In such systems, the two-point function shows a 
power- law decay just as in our problem [22|. The effect 
of such long range potentials on the critical exponents 
have been studied earlier [13, US [l^-HB] with the general 
conclusion that the long range part of the interaction do 
influence the critical behavior of the system More 
precisely, if we have an attractive pair potential in d di- 
mensions of the form —4'{r) ^ ^^^^^ a > 0, then 
critical exponents take their short-range values for all 
(T > 2 — rjsr- For cr < 2 — rj^r, two kinds of behavior exist. 
For < a < d/2,the exponents take their mean-field val- 
ues and for d/2 < a < 2, the exponents depend on the 
value of a (See cite aizenman and references therein). 

Though there is a well established connection between 
the lattice percolation problem and the Ising model [111 , 
there is no similar result connecting the continuum perco- 
lation problem to any simple Hamiltonian system. How- 
ever, simulation studies suggest that the continuum per- 
colation problem (where percolating units are bounded 
in size) and the lattice per colation problem belong to the 
same universality class 15]. In our problem, the behavior 
is similar to that of an Ising system with long range po- 
tential of the form —(t>{r) ^ where a > 2. Then the 

results from 27[ shows that a deviation from the stan- 
dard critical behavior is expected when a < 3 — rjsr and 
the critical exponents will take their short-range values 
for a > 3 — rjsr- Also mean- field behavior is expected 
when a < 2. However for this range of a, the entire plane 
is covered for all non-zero number densities and hence 
there is no phase transition. 

In the next two sections, we investigate for the de- 
pendence of exponents on the power-law tail of the dis- 
tribution of the radii of the discs. First we develop an 
approximate RG method. Then we carry out simulation 
studies which show that the correlation length exponent 
V takes its short range value for 
depends upon a for a < 3 — rjsr ■ 



> 3 — rjsr, while it 



IV. CRITICAL BEHAVIOR NEAR THE 
PERCOLATION THRESHOLD 

In this section, we propose an approximate RG method 
to analyze the behavior of continuum percolation mod- 
els near the percolation threshold, when the percolating 
units have a distribution of sizes. We assume that we can 
replace discs of one size having a number density n with 
discs of another size and number density n', provided 
the correlation length remains the same. Application of 
a similar idea in disc percolation problem with only two 
sizes of discs may be found in Q . 

We will illustrate the method by considering a prob- 
lem in which the radii of discs take only two possible 
values, say i?i and i?2- Let their areal densities be Ai 
and A2 respectively, and assume that both Ai and A2 are 
below A* , the critical threshold for the percolation prob- 
lem with only single sized discs present { A* 1.128085 
p^V Also let ^1 represent the correlation length when 
only discs of size Ri are present in the system and ^2 
represent that when only discs of size R2 are present. In- 
variance of the two point function under length rescaling 
requires that the expression for the correlation length ^ 
is of the form ^ — Eg (A), where the function g{A) de- 
termines how the correlation length depends on the areal 
density A and is independent of the radius R. Let A2 is 
the areal density of the discs of size R2 which will give 
the same correlation length as the discs of size Ri. i.e. 



(9) 



or 



a (^i) = 6U2 



Rig{Ai)^R2g[A 



(10) 



Given the form of the function g{A), we can invert the 
above equation to find A2. Formally, 



A2 = 9-' 



(11) 



So the problem is reduced to one in which only discs of 
size i?2 are present, whose net areal density is now given 
by. 



A'^A^+ An 



(12) 



System percolates when A'2 = A*. Now, when areal den- 
sity A is close to A*, we have 



g{A) ^ K {A* - Ay" 



(13) 



where K is some constant independent of A and v is 
the correlation-length exponent in the usual percolation 
problem. Using this in Eq. [TI] we get 



A2 = A*- {A* - Ai) {R2/R1) 



(14) 



Therefore, for a given value of Ax < A* , the areal density 
of discs of radius i?2, so that the system becomes critical 
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Figure 2: Variation of £m with m for different values of po 
showing sub critical and supercritical regimes. We have used 
a = 3 and A = 2. 



is given by, 



A. 



{A* -A,)iR2/Ri) 



(15) 



So the total areal density at the percolation threshold is, 

+ ^2 = + (A* - Ai) {R2/Rif''' 
= Ai{\ - x) + A*x 



where x — {R2/R1) 



Without loss of generality we 
may assume i?2 > R\- Then a; > 1 and we can see 
from the above expression that the percolation threshold 
Ai + ^2 > ^* , a result well known from both theoretical 
studies [1^ and simulation studies |l3| . 

Now in our problem assume that areal density of discs 
of type do not exceed A*. Renormalizing discs up to 
type m in our problem gives the equation for the effective 
areal density of the m-th type discs as 



A^ 



A* 



{A* ~K^_^)h'l^ 



(16) 



where m > 1, Aq = po and p,„ = noTrA^^""^"* denote 
the areal density of discs of radius A™. Here no is the 
number density of discs of radius i?o (or of type 0), which 
for convenience we have set equal to unity. If we denote 
A* — A'j^ by Em which is the distance from the criticality 
after m-th step of the renormalization, then the above 
expression becomes 



Pn 



(17) 



The equation describes the flow near the critical point 
when we start with a value of pQ, the areal density of the 
first type of discs. Here gives the effective distance 
from criticality of the m-th order discs in the system, in 
which now only m-th and higher order discs are present. 
Now for given values of the parameters a and A, we can 
evaluate Sm in Eq- [13 using a computer program and 



plot Em verses m. Depending upon the value of po, we 
get three different behaviors. For value of po below the 
critical value denoted by pg , e„i will go to A* asymptoti- 
cally (System is sub critical) and when it is above Pq, Sm 
will go to —00 asymptotically (System is super critical). 
As Po — > Po , we get the critical behavior characterized by 
Em tending to the RG fixed point asymptotically. Typ- 
ical result using Eq [T7] with A = 2 and a = 3 is shown 
in Fig [21 We can see that as we tune po, the system 
approaches criticality, staying closer to the Sm — line 
longer and longer. Critical behavior here can be char- 
acterized by the value of m at which the curve deviates 
from the approach to Sm = line. To understand how 
the correlation length diverges as we approach critical- 
ity, we assume that we can replace the sub critical sys- 
tem with a system where only discs of type m' is present 
and has a fixed areal density below A*, where m' is the 
value of m at which Em shows a substantial increase - say 
Em becomes A*/2. For continuum percolation problem 
with single sized discs, the correlation length ^ = Rg{A), 
where g{A) is a function with no explicit dependence on 
radius R. Therefore, correlation length in our problem. 



^ a A" 



(18) 



We can write the recurrence relation Eq. (|17p in terms 
of the areal density p„ as 



(19) 



n=0 



But p„ = poA"(2-a)^ Therefore, 



1 - M^-a-l/f)] 



(20) 



For large values of m, the last term in the above equation 
involving J^^n{2-a-i/u) |-jg neglected. Then, 



aI^] 



A* 



Po 



1 - A(2-a-l/!^) 



Therefore, 



a[^] 



Po 



1 - A(2-a-l/!^) 



(21) 



(22) 



For a given value of po < A*, the order m' at which Em 
is increased substantially, say to a value A* /2 is given by 



^' = [logA {A*/2) - logA {p*o ~ po) 
-KlogA(l-A(2-'^-i/'-)^]j,+ l 



(23) 



So for Po close to pp and large values of a, 

m' ^ logA (pS - Po)"'' . (24) 
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Figure 3: Variation of n* with a for two different values of A 
(See Eq. 128(1 . The horizontal line corresponds to the threshold 
for the single sized discs case. 



so that 



(25) 



Thus we find that the correlation length exponent v is 
independent of the parameters a and A of the distribu- 
tion. From Eq. [2H we can also obtain the percolation 
threshold pp as a function of the parameters a and A. In 
Eq. [22] left hand side is positive definite. So for values 
of pq for which j^^j^^J'-a-i/i,) < A*, we will have Em > 
for large values of m. Similarly for values of po for which 
^_^(2^°„_iy„) > A*, we will have < for large values 
of m. Hence the critical areal density must be given 
by 



A* = 



Po 



[1 _ A(2~a-l/i')] 



(26) 



Therefore, 



p*^A* 



Or in terms of the total number density, the percolation 
threshold n* is given by. 



(28) 



where Uc = A*/^, the critical number density for perco- 
lation with single sized discs of unit radius. Note that 
this approximate result does not give the correct limit, 
n* — >■ as a — >■ 2. The RG scheme depends on the ap- 
proximation that the effect of size i?i of areal density Ai 
is the same as that of discs of radius i?2 of density A2 , as 
in Eq. [Q] This is apparently good only for a > 3 — rjsr- 
Fig [3] shows the variation of the critical threshold with a 
for two different values of A using Eq. [^S] As one would 
expect, for large values of a, n* tends to ric- 



V. SIMULATION RESULTS 

We determine the exponent v and the percolation 
threshold n* by simulating the continuum percolation 
system in 2D, with discs having a power law distribution 
for their radii. We assume a continuous distribution for 
the radii where given a disc, the probability that it has a 
radius between R and R + dR is equal to aR^^°-'^^'> where 
a > 2. For a < 2 the entire plane is covered for arbitrarily 
low densities of the discs. We use cyclic boundary con- 
ditions and consider the system as percolating whenever 
it has a path through the discs from the left to the right 
boundary. We drop discs one at a time on to a region of 
a plane of size L x L, each time checking whether the sys- 
tem has formed a spanning cluster or not. Thus number 
density is increased in steps of So after dropping 

the n — th disc, the number density is n/L^. Now associ- 
ated with each number density we have a counter say /„ 
which is initialized to in the beginning. If the system is 
found to span after dropping the n'-th disc, then all coun- 
ters for n > n' is incremented by one. After a spanning 
cluster is formed, we stop. By this way we can determine 
the spanning probability n{n, L) — fn/N where N is the 
number of realizations sampled. The number of realiza- 
tions sampled varies from a maximum of 2.75 x 10^ for 
a = 2.05 and L = 90 to a minimum of 4000 for a = 10.0 
and L = 1020. This method of dropping basic percolat- 
ing units one by one until the spanning cluster is formed 
has been used before [2^ in the context of stick percola- 
tion which was based on the algorithm developed in (30| , 
and allows us to study relatively large system sizes with 
large number of realizations within reasonable time. 

The probability that there is at least a single disc which 
span the system of size L at number density n is 1 — 

^exp~"^° /L°. It is easy to see that to leading order 

in n, this 'long range 'part of the spanning probability 

n2°' ~ 



n{n,L)ir is 



So one can write a scaling form for 



L"-2 

the spanning probability, 
(27) n{n, L) = i7(n, L)ir -f (1 - 7T(n, L)j,)(/.((n* 



(29) 

Therefore we can define the 'short range 'part of the span- 
ning probability n'{n,L) = {n{n,L) - n{n,L)ir)/{l - 
n{n, L)ir) where the leading long range part is sub- 
tracted out. Therefore, we have 



n'{n,L) = (f){{n* -71)1^/") 



(30) 



and the scaling relations, (See for eg. [31j ) 

A(L) a L-^/" (31) 



a A 



(32) 



where n*j:j:(L) is a suitable defined effective percolation 
threshold for the system of size L, and A is the width 
of the percolation transition obtained from the spanning 
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Figure 4: Plot of effective percolation threshold nlff against 
A for a — 2.25 and a — 3.25. The best straight line fit is 
obtained with the last four data points. 





Log(L) 



Log(L) 




Figure 6: Variation of l/u with a. The horizontal line corre- 
sponds to the standard 2D percolation value l/v = 3/4. 




Figure 5: Log-Log plot of A Vs L for a — 2.25 and a = 4.0 
along with lines of slope —.47 and —.75. 



Figure 7: Variation of percolation threshold n* with a. The 
horizontal line corresponds to the threshold for the single sized 
discs case. (Inset) Asymtotic approach of n* to the single 
sized discs value Uc — .3591 along with a straight line of slope 
-1. 



probability curves n'{n, L). Note that Eqs. [3l1and[32lare 
applicable with any consistent definition of the effective 
percolation threshold and width A [3l|. A good way to 
obtain n*j:j: and A is to fit the sigmoidal shaped curves 
of the spanning probability n'{n,L) with the function 
1 / 2[1 + er f [{n-nl J ^{L))/A{L)]] (see [H), which defines 
the effective percolation threshold "(^e// the number 
density at which the spanning probability is 1/2. We 



determined 



and A for each value of a and L and 



determined l/v and n* for different values of a using Eqs. 
1311 and [32] respectively. Typical examples are shown in 
fig. Eland Fig. O 

At first, we determined the percolation threshold and 
the exponent for a system of single sized discs of unit ra- 
dius. We obtained n* = .3589(±.0001) (or areal density 
w 1.12752) and = .758(±.018) in very good agree- 
ment with the known value for the threshold TS] and the 



conjectured value of 1/z^ = 3/4 for the exponent. Values 
of l/v obtained for various values of a are shown in fig|6l 
We scan the low a regime more closely for any variation 
from the standard answer. We can see that the estimates 
for 1/v are very much in line with the standard perco- 
lation value for a > 3 — rjar while it varies with a for 
a < 3 — rjsr- Fig. [7] shows the variation of the percolation 
threshold n* with a. As expected, with increasing a, the 
percolation threshold increases and tends to the single 
sized disc value as a — )■ oo, and as a — >■ 2, the thresh- 
old tends to zero. Tha data also shows that n* coverges 
to the threshold for the single sized disc value as 1/a. 
Values of the threshold for some values of a are given in 
Table H 

Finally as a check, we plot the spanning probability 
n'{n,L) (see Eq. ^ against (n - n*)L^^'' to be sure 
that a good scaling collapse is obtained. We show two 
such plots for a = 2.50 and a = 4 in fig. [S] We can 



7 



a 


n* 


r]* = n*na/{a — 2) 


= 1 — exp 


2.05 


0.0380(6) 


4.90(7) 


0.993(1) 


2.25 


0.0693(1) 


1.959(3) 


0.8591(5) 


2.50 


0.09745(11) 


1.5307(17) 


0.7836(4) 


3.50 


0.16679(8) 


1.2226(6) 


0.70555(17) 


4.00 


0.18916(3) 


1.1885(2) 


0.69543(6) 


tJ.UU 




1 1 '^Q7(A\ 


n P,RP,d'^^^ 


6.00 


0.24340(5) 


1.1470(2) 


0.68241(8) 


7.00 


0.2593(2) 


1.1406(7) 


0.6804(2) 


8.00 


0.27140(7) 


1.1368(3) 


0.67917(9) 


9.00 


0.28098(9) 


1.1349(4) 


0.67856(12) 



Table I: Percolation threshold n* for a few values of a along 
with corresponding critical areal density rj* and the critical 
covered area fraction 6* . 




a-4.00 


1 7 




\ L = 200 




/ L = 300 




/ L = 450 




L = 850 




L- 1020 



0.03 0.06 0.09 0.12 



n =.09745 
1/v=.65 




VI. CONCLUDING REMARKS 

In this paper, we discuss the effect of a power-law 
distribution of the radii on the critical behavior of a 
disc percolation system. If the distribution of radii is 
bounded, then one would expect the critical exponents 
to be unchanged. However, if the distribution of radii 
has a power-law tail, we show that this strongly influ- 
ence the nature of the phase transition. The whole of 
the low-density non-percolating phase has power-law de- 
cay of correlations. And this occurs for any value of the 
power a, howsoever large. The critical exponents depend 
on the value of a for a < 3—r]sr and take their short-range 
values for a > 3 — rysr- We also propose an approximate 
RG scheme to analyse such systems. Using this, we com- 
pute the correlation-length exponent and the percolation 
threshold. The approximate RG scheme is good only for 
a > 3 — rjsr- Monte-Carlo simulation results for the per- 
colation thresholds and the correlation-length exponent 
are presented. 

We can easily extend the discussion to higher dimen- 
sions, or other shapes of objects. It is easy to see that the 
power law correlations will exist in corresponding prob- 
lems in higher dimensions as well. 



(n - n')L"" 
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Figure 8: Variation of il(n, L) with n (top row) and the scal- 
ing collapse (bottom row) for a — 2.50 (on left) and a = 4.00 
(on right). 

see that a very good collapse is obtained. Similar good 
collapse is obtained for other values of a as well. 
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